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Abstract 



A search for 77 decays of a Higgs boson is performed in the data sample collected at LEP 
with the ALEPH detector between 1991 and 1999. This corresponds to an integrated luminosity 
of 672 pb~ x at centre-of-mass energies ranging from 88 to 202 GeV. The search is based on 
topologies arising from a Higgs boson produced in association with a fermion pair via the Higgs- 
strahlung process e + e" — > Hff, with ff = uu, e + e~, fi + fi~, t + t~ or qq. Twenty-two events are 
selected in the data, while 28 events are expected from standard model processes. An upper 
limit is derived, as a function of the Higgs boson mass, on the product of the e + e~ — > Hff cross 
section and the H — > 77 branching fraction. In particular, a fermiophobic Higgs boson produced 
with the standard model cross section is excluded at 95% confidence level for all masses below 
100.7 GeV/c 2 . 
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1 Introduction 



In general, neutral Higgs bosons do not couple directly to massless photons. For instance 
the standard model Higgs boson couples to photons only through loops of charged particles, 
i.e., W's, quarks and leptons, and the branching ratio into 77 is small (~ fl| for 

hih ~ 90GeV/c 2 ). However the 77 branching fraction of Higgs bosons can be increased with 
respect to the standard model prediction in each of the following four configurations. 

• The direct couplings to fermions are suppressed, as is the case for models with at least two 
Higgs multiplets [fj, of which one couples only to fermions and the others only to gauge 
bosons. The physical states with couplings only to gauge bosons are called fermiophobic 
Higgs bosons. 

• The direct couplings to gauge bosons are enhanced with anomalous couplings 0. These 
couplings are described in the most general formulation with four effective six-dimensional 
operators with strength fj/A 2 , where A is the scale of the new underlying interaction. 

• Couplings to both fermions and bosons are modified, as is the case in the minimal 
super symmetric extension of the standard model (MSSM). 

• Additional light, charged particles enter the loops that couple Higgs bosons and photons, 
as is again the case in the MSSM (charginos, squarks, sleptons, charged Higgs bosons). 

With some particular choices of parameters, the branching ratio into 77 may be enhanced in 
the MSSM, and can reach a value close to 100% in models with fermiophobia or with anomalous 
gauge couplings. It is therefore possible that a Higgs boson has escaped the standard search 
for the Higgs-strahlung process e + e~ — ■> Hff with H — > bb [|]]. In this letter, a complementary 
search for the Higgs-strahlung process with H —>■ 77 is described. 

The analysis addresses all topologies arising from the e + e~ — > HZ^ process, as characterized 
by the charged track multiplicity of the final state: (i) acoplanar photons with missing energy 
and no charged particles for B.uu; (ii) photon pairs with exactly two charged particles identified 
as leptons for H£ + £~; (Hi) photon pairs accompanied with two thin, low multiplicity jets for 
Hr + r~, from two to four charged particles; and (iv) photon pairs with a hadronic system for 
Hqq with at least five charged particles. 

The analysis is performed with the data collected with the ALEPH detector from 1991 to 
1999 including the Z peak data collected during the LEP 2 period. This sample corresponds to 
an integrated luminosity of 672pb _1 at centre-of-mass energies ranging from 88 to 202 GeV. 
Details are given in Table [3]. 

After a short description of the detector properties relevant for this search, the common 
preselection based on photon identification is reviewed in Section 3. The global search strategy 
is developed in Section 4. The systematic uncertainties affecting the selection efficiency are 
discussed in Section 5 and the results are given in Section 6. 
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Table 1: Integrated luminosity for the Z peak data (88-94 GeV) and high energy data (130-202 GeV). The 
energies are rounded to the closest integer value. 



(GeV) 


88 


89 


90 


91 


92 


93 


94 


130 


136 


L (pb- 1 ) 


0.7 


17.9 


0.8 


124.8 


0.8 


19.1 


0.8 


6.2 


6.4 


(GeV) 


161 


170 


172 


183 


189 


192 


196 


200 


202 


L (pb- 1 ) 


11.1 


1.1 


9.5 


59.2 


177.1 


28.9 


79.8 


86.3 


42.0 



2 The ALEPH detector 

The ALEPH detector and its performance are described in Refs. f|||. The tracking detectors, 
composed of the silicon vertex detector surrounded by the inner tracking chamber and the time 
projection chamber (TPC), provide efficient reconstruction of charged particles in the angular 
range |cos#| < 0.96. A charged particle track is called a good track if it is reconstructed with 
a least four hits in the TPC and if it originated from within a cylinder of length 20 cm and 
radius 2 cm, coaxial with the beam and centred at the interaction point. A 1.5 T axial magnetic 
field delivered by a super-conducting solenoidal coil allows a charged particle l/p_i_ resolution 
of (6 x 10~ 4 ©5 x 10- 3 /p ± )(GeV/c)- 1 to be achieved. 

The electromagnetic calorimeter (ECAL) is a lead/ wire-plane sampling calorimeter covering 
the angular range |cos#| < 0.98. Anode wire signals provide a measurement of the arrival time 
of the particles relative to the beam crossing with a resolution better than 15 ns. Cathode pads 
associated with each wire layer are connected to form projective towers of approximately 0.9° 
by 0.9° which are read out in three segments in depth. The impact parameter of the photon 
with respect to the interaction point is estimated from the barycentre of the electromagnetic 
shower in each segment with a resolution of about 6 cm. A photon candidate is identified 
using a topological search |G| for energy deposits in neighbouring electromagnetic calorimeter 
towers isolated from the extrapolation of any charged particle track to the ECAL. Any photon 
candidate close to a boundary between ECAL modules or pointing towards an uninstrumented 
region of the TPC is not considered in the analysis. The energy calibration of the ECAL is 
obtained from Bhabha events, radiative returns to the Z resonance, e + e~ — > 77 and 77 — > e + e~ 
events. The energy resolution for photons is <5E/E = 0.25/^E/GeV + 0.009 [[]]. 

The luminosity monitors (LCAL and SICAL) extend the calorimetric coverage down to 
small polar angles. The iron return yoke is instrumented with streamer tubes and acts as a 
hadron calorimeter (HCAL), covering polar angles down to llOmrad. Surrounding the HCAL 
are two additional double layers of streamer tubes called muon chambers. 

The measurements of the tracking detectors and the calorimeters are combined into objects 
classified as charged particles, photons and neutral hadrons using the energy flow algorithm 
described in Ref. ||. All objects are used to compute the total visible energy E vis and the 
missing energy E miss with a resolution of (0.6-y/E vis /GeV + 0.6) GeV. 

Electron identification is based on the matching between the measured momentum in 
the tracking system and the energy in the ECAL, the shower profile in the ECAL and the 
measurement of the specific ionisation energy loss in the TPC. Muons are identified by their 
characteristic hit patterns in the HCAL and in the muon chambers. 
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3 Event preselection 



Signal events are characterized by two isolated, energetic photons, well contained in the 
apparatus, and in time with the beam crossing. 

Photon isolation is ensured by requiring that the total charged energy in a cone of half- 
angle 14° around the photon direction (E^o ) be smaller than 2 GeV, and that the invariant 
mass between the photon and any charged particle (m 7jC h) be in excess of 1 GeV/c 2 . Pairs of 
photons from a 7r° decay are rejected by the requirement that their invariant mass be larger 
than 1 GeV/c 2 . 

Photon centrality is enforced by requiring that the polar angles satisfy |cos# 7 | < 0.9. Since 
the Higgs boson is expected to be produced with a nearly uniform cos 9 distribution, the sum 
of the two cosines (in absolute value) must be smaller than 1.4, and the cosines of the Higgs 
boson production and the thrust axis polar angles, # 77 and thrust must be within ±0.95. Events 
with energy at low polar angles are further vetoed by rejecting events with more than 2 GeV 
within 14° of the beam axis (E\la m ). 

If, in a selected event, more than two photons satisfy the above criteria, only the most 
energetic two photons are considered as originating from a Higgs boson candidate. 

Finally, only events in time with the beam crossing are kept. For events with at least 
two good tracks, a timing to better than 1 ns is ensured by the good track definition. Events 
with only one good track are rejected. The production time of events with no good tracks is 
determined as the energy-weighted average of the times reconstructed in all ECAL modules, 
to, and is required to be in agreement with the beam crossing time within ±40 ns. For these 
events, the impact parameter of each identified photon is also required to be less than 25 cm. 

The reconstructed particles of events with at least two good tracks are then forced to form 
four "jets" with the Durham jet clustering algorithm J/J. The consistency of these events with 
a four-body final-state hypothesis is verified by requiring that each of the four jet energies Ef sc , 
rescaled to satisfy energy-momentum conservation under the assumption that the jet velocities 
are perfectly measured, be positive. However, in order to make the reconstructed Higgs boson 
mass resolution independent of the final state topology, the measured photon energies were 
used instead of the rescaled energies. 

In the following, the neutral electromagnetic energy in a given jet is computed with all neutral 
objects in the electromagnetic calorimeter, and with the objects in the hadron calorimeter found 
behind uninstrumented regions of the electromagnetic calorimeter. An electromagnetic jet is 
defined as a jet with more than 80% of electromagnetic energy. In events reconstructed with 
four jets, i.e., with at least two good tracks, the two jets with the smallest electromagnetic 
energy fraction are called fermionic jets. 

At this level, a signal efficiency of 20 to 65% is achieved at all centre-of-mass energies and 
for any Higgs boson mass above lGeV/c 2 , as estimated with many e + e~ — * Hff simulated 
event samples produced with the HZHA generator || and processed through the whole detector 
simulation and event reconstruction chain. 
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4 Event selection 



Events passing the above preselection are further classified in the four signal topologies 
according to their good track multiplicity (n c h). 



Events with no good tracks are classified as Hz/i/ candidate events. The main standard 
model background sources to this final state are (i) e + e~ — > 1^7(7), simulated with the 
KORALZ package 0; and (%%) e + e~ — > 77(7), simulated with the GGG generator |l(J. The 
latter generator does not contain QED contributions of order a 4 and above, but they 



were estimated in Ref. [IT| to be small enough to be considered negligible in the present 
analysis. 

Events with two good tracks, both positively identified as electrons or muons, are classified 
as H£ + £~ (with I = e or //) candidate events. The programs BHWIDE [12| and UNIBAB [13 



on the one hand, and KORALZ on the other, are employed to simulate the main background 



processes, i.e., e T e 



e + e ry(>y) and 7(7), respectively. 



Other two good track events, and events with up to four good tracks are classified as 
Ht + t~ candidate events. Here again, the main background process, e + e~ — > t + t~7(7) 
is simulated with KORALZ. 

Finally, events with at least five good tracks are classified as Hqq candidate events. The 
e + e~ 



• 0*1(7) process is simulated with JETSET [14| for data taken at the Z resonance, 
and with PYTHIA [15j and HERWIG ]TB| at higher energies. The four-fermion processes 



WW, ZZ, Zee and Wez/ are simulated using PYTHIA |L5| 



The selection criteria designed to reduce the contribution of the background processes are 
summarized in Table [2], and only a brief account is given here. 

For each of the topologies, the photons of the e + e~ — > 1*17(7) background process are emitted 
by the incoming and the outgoing charged fermions. They are therefore preferentially produced 
either along the beam axis or along one of the outgoing charged fermion directions. An efficient 
rejection is achieved by tightening the photon isolation requirements with respect to these 
directions, by means of the variables introduced in Section 3 and a number of other relevant 
variables: (i) photon energies E 7 and the energy recoiling against the photon pair EJ| C011 , (ii) 
invariant masses (m 7jC h, m 7i f) and transverse momentum pT of individual charged particles 
and fermionic jets with respect to the photons, (in) transverse momentum of the photon with 
respect to the closest fermionic jet p^ and transverse momentum of the fermionic jets with 
respect to the thrust axis p^ f rust , (iv) the value y 3 4 of the Durham y cut transition value between 
three and four jets and the invariant mass of the two fermionic jets m f j. For the Z peak data, 
exactly two electomagnetic jets Nj^s are required in order to reject events in which both photons 
are inside the same jet. 

Further cuts are made on the cosine of the decay angle of the photons (fermions) in the 
rest frame of the 77 (ff) system cos^ 77 y (cos6^ cy ). These decay angles are expected to have 
a flat distribution for signal events, while they are strongly peaked towards small angles for 
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Table 2: Overview of all selection cuts 



Z peak 



High energy 



E 7i > 3GeV; E™ ax > O.ly/s 
Eg"* 1 > 10 GeV; E visiblc > 0.6^ 



0.2E 



beam 



< E 7 < 0.75E beam 



|m rec — m z | < 15 GeV/ c 2 



No good track:Yivv 



I cos 9. 



77 



m 



E miss > 15 GeV 
E vis - E 77 < 10 GeV 
-, 7 + m rec — y/s\ > 2 GeV/ c 2 
Icosgg l < 0.95 



< 0.9 



E max > 



Two good tracks:H£ + £ /Hr + r 



pl ch > 2 GeV/c 
m 7i f> 10 GeV/ c 2 



Ef sc > GeV 
m ff > 1 GeV/c 2 
m 7i ch > 5 GeV/c 2 



e+e'/fi+fi' : E£+t 



E 7 > 10 GeV 
Icosflffi-I < 0.95 



e/z/eh//di/hh : Hr + r 



E 7 > 20 GeV 



— m 



ff 



< 20 GeV/c A 



pl ch >10GeV/c 



Three or four good tracks:B.T + r 



Er c > GeV 



pthmst > 2 GeV/c; m. 
Ng = 2; E 



7if > 10 GeV/ c 2 
< 15 GeV 



cos^ 7 y | < 0.95; |cos^ cy | 



< 0.95 



pthmst > 4GeV/c 



Five or more good tracks:Hqq 



Ef sc > GeV 
|cos# 7 | < 0.8 

?/34> 0.001 



„thrust 
P-Lf 



> 2GeV/c;m 7if > 10 GeV/ c 2 
Nj s c = 2; E miss <15GeV 
cos^ 7 y | < 0.95; |cos^ cy | < 0.95 



pthmst > 4GeV / c 
pf > 0.05^ 



i_ 7 

^min 



4— jets 



> 350° 



the background. For a high mass Higgs at high energies the signal events are expected to 
be somewhat spherical in nature. Non-spherical events are rejected by a cut on the sum of 
the four minimum inter-jet angles 9™^ cts . At high energies the mass recoiling against the 
photonic system m rec is required to be consistent with the Z mass. For the Hrr channel at 
Z peak energies the mass recoiling against the photonic system is required to be consistent 
with the rescaled mass of the fermionic system mff™. Finally, in the B.uu topology, some boost 
is ensured by requiring that the total produced mass does not saturate the available energy: 
|m 77 + m rec — y/s\ > 2 GeV. 

The typical efficiencies and expected standard model background for each topology, together 
with cross-channel contamination, are shown in Table |3]. 
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Table 3: Efficiencies for the different Z decay channels for a 110GeV/c 2 Higgs boson mass at *Js = 202 GeV 
and numbers of expected background and observed events at the Z peak (Z) and at high energies (HE). 



n c h 




£m+e- 




CHqq 


N^f (Z) 


N obs .(Z) 


N exp k ' (HE) 


N obs . (HE) 





47.0% 


0.3% 


0.0 


0.0 


0.9 


1 


3.6 





2 


0.0 


36.7% 


22.4% 


0.0 


2.3 


4 


2.9 


2 


3-4 


0.0 


2.0% 


14.9% 


0.0 


0.7 





0.4 


2 


>5 


0.0 


0.0 


5.2% 


37.3% 


9.3 


7 


7.8 


6 


Total 


47.0% 


39.0% 


42.5% 


37.3% 


13.2 


12 


14.7 
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5 Systematic uncertainties 



Because no background subtraction is performed to derive the final result, the uncertainty on 
the background evaluation has not been estimated. The uncertainty on the photon selection 
efficiency summarized in Table ^ receives a large contribution from the photon identification 
efficiency. The uncertainties due to photon energy calibration or photon angle resolution affect 
only the 77 invariant mass resolution. 

To estimate the uncertainty on the photon selection efficiency, the total cross section of the 
e + e~ — > 77(7) process was measured and compared to its prediction. Two-photon-candidate 
events are selected as final states with only two identified photons (with the same timing and 
pointing constraints as those defined in section 3) with a polar angle such that |cos# 7 | < 0.9 
(0.95), and an opening angle satisfying cosa 77 < —0.999 (—0.9999), for the Z peak data (at 
high energy). The e + e~ — > e + e~77 background is rejected by requiring no charged particles in 
the event and not more than 2 GeV of energy around the beam direction (E^ 111 < 2 GeV). 
The ^z/77 background is reduced to a negligible amount by requiring the invariant mass of the 
two-photon system to be greater than 0.75-y/s. The result is displayed in Fig. [I]. From this 
measurement, a conservative relative systematic uncertainty is estimated to be 3%, including 
a 1% theoretical error. The systematic uncertainty on the photon detection efficiency includes 
the effects of the cuts on photons used at preselection level as well as that of the integrated 
luminosity determination, which enters the e + e~ — > 77(7) cross section measurement. 

The effect of the photon energy calibration and angular resolution on the efficiency is 
estimated with Bhabha events. The relative difference between the values of the electron energy 
determined independently by the tracking system and by the calorimetry mainly originates from 
electron Bremsstrahlung. This difference is measured to be ~ 1.5%, and it agrees within 0.2% 
with that expected from the simulation. This ±0.2% uncertainty is conservatively assigned 
to the photon energy calibration. It leads to a shift of the measured Higgs boson mass 
by ±200 MeV/ c 2 , and to a negligible increase of the mass resolution. The uncertainty on 
the reconstructed photon direction is estimated by comparing the directions determined with 
the tracking and the calorimetry. It degrades the measured Higgs boson mass resolution by 
approximately 20 MeV/c 2 with a negligible loss of efficiency. These systematic uncertainties 
are included in the final result by increasing the Higgs boson mass resolution by 200 MeV/ c 2 . 

As the lepton identification is used only to classify events, the signal efficiency decreases by 
less than 0.1% relative when the lepton identification efficiency is modified by 10%. 
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Figure 1: Ratio of the measured cross section of the e + e~ — > 77(7) process at all centre-of-mass energies to 
the expected standard model cross section. The solid line is the best fit to the data and the two dotted lines 
represent the statistical error from the fit. 



The remaining selection criteria are expected to be largely insensitive to the details of the 
simulation of the hadronic system. A quantitative estimate of the size of possible discrepancies 
is performed with the aid of an event reweighting technique. For each selection variable, bin-by- 
bin correction factors are calculated as the ratio of data to Monte Carlo expectation, evaluated 
at the preselection level with the cut on m 7jC h removed. The Monte Carlo signal distribution 
of the selection variable, obtained when all cuts are applied, is then re-weighted with these 
correction factors to obtain a new estimate of the efficiency The difference in efficiencies 
estimated with this technique are, for almost all variables, at the level of a few parts per mil 
and are given in Table |j. The largest effect, of 1.2%, comes from the cut on the isolation of the 
photon m 7jC h. 

Finally, the model dependence of this analysis is estimated with the anomalous coupling 
model mentioned in Section 1. In that context, all the anomalous coupling parameters fj/A 2 are 
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Table 4: Systematic uncertainties on signal efficiency 



Sources 


Relative uncertainty in % 


Photon selection efficiency 


3.0 


7 energy calibration 


0.5 


7 angular resolution 


0.1 


Total energy calibration (Z peak data only) 


0.5 


Lepton identification 


0.1 


Photon isolation 


1.2 


COS ^thrust 


0.1 


.^thrust 
P-U 


U.Z 


m 7,f 


0.3 


Pj_^ (High energy data only) 


0.4 


Y34 


0.5 


^T-\cts (High energy data only) 


0.3 


Model dependence 


4.0 


Total in quadrature 


5.2 



varied independently of each other within ±100 TeV -2 using the HZHA generator. The influence 
of the Higgs boson energy and angular distributions leads to a relative uncertainty on the signal 
efficiency of at most 4%. 

The decay width of the Higgs boson is in general negligible with respect to the detector 
resolution (below 1 GeV for |fj/A 2 | < 100 TeV™ 2 ). In anomalous coupling models with |fj/A 2 | 
above 500 TeV~ 2 the decay width of a heavy Higgs boson becomes larger than the detector 
invariant mass resolution (r ~ 3 GeV) and the relative variation of the signal efficiency becomes 
larger than the previously estimated systematic uncertainties. The final results presented in 
the next section are therefore valid for models in which the width of the Higgs boson is less 
than a few GeV. 



6 Results 

The signal efficiencies are displayed in Fig. |a as a function of the Higgs boson mass hypothesis. 
Twenty-two events are selected in the data compared with an expectation of 27.9 ±1.6 events 
from all standard model background processes. The contributions of the various channels are 
given in Table The diphoton invariant mass distribution of these selected events is shown for 
both data and Monte Carlo expectation in Fig. ^b. The low mass background is dominated by 
events where both photons originated from the same jet, whereas the high mass background is 
dominated by events where the two photons originated from different jets. The sum of these 
two contributions has a minimum at around 8GeV/c 2 . No evidence of a resonance decaying 
to 77 is observed and a 95% confidence level upper limit on the number of signal events at 



a given diphoton invariant mass is derived following the method described in Ref. ||17|| . For 
each LEP energy, the total efficiency and the mass resolution are parametrized as a function of 
the 77 invariant mass. The diphoton invariant mass resolution varies linearly from 1 GeV/c 2 
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Figure 2: (a) Hff efficiency at the Z peak and high energy (yfs — 202 GeV). (b) Diphoton invariant mass 
distribution for data (dots with error bar), expected sources of background (solid histogram) and a 110 GeV/ c 2 
Higgs boson signal at -^s = 202 GeV with arbitrary normalization (dashed histogram) . 



to 3.5 GeV/ c 2 over the whole Higgs boson mass range. The systematic uncertainties are 
conservatively taken into account by scaling down the efficiency by 5.2% and increasing the 
mass resolution as described in the previous section. 

The present analysis does not apply for Higgs boson masses below 1 GeV/ c 2 . To extend the 
search to smaller masses, the direct measurement of the Z invisible width with single photon 



counting described in Ref. ITS! was used. For this measurement, the single photon candidate 
events were selected as final states with only one cluster in the electromagnetic calorimeter. 
The Huu events with hih < 1 GeV/ c 2 would therefore have been selected with a good efficiency, 
ranging from 13 to 45%, thus leading to a sensitivity similar to that of the present analysis. 

Assuming a Higgs boson production cross section with the same \fs, me dependencies as 
in the standard model, the 95% confidence level upper limit on the product branching ratio 
B(H — > 77)<j(e + e - — > Hff)/cr SM (e + e _ — > Hff) is derived and shown in Fig. |3|. 

For the case of a Higgs boson produced at the standard model rate, the best upper limit on 
the branching ratio (4.7x 10~ 3 at 95% confidence level) is obtained for Higgs boson masses below 
20 GeV/ c 2 . A Higgs boson decaying exclusively to two photons is ruled out up to 109 GeV/ c 2 
at 95% confidence level. 



A fermiophobic Higgs boson |19[ with no tree-level coupling to fermions is excluded at 95% 
confidence level for any mass up to 100.7 GeV/ c 2 . 

The present analysis extends the reach of similar analyses performed by other LEP 
collaborations [E0TE21. 
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Figure 3: Measured (full curve) and expected (dash-dotted curve) 95% confidence level upper limit 
on B(H — > 77)<r(e + e - — > Hff)/cr SM (e + e~ — > Hff). The dashed curve is the predicted branching ratio for a 
fermiophobic Higgs boson in the limit of B(H — > ff) = 0. 

7 Conclusion 



With a data sample of 672 pb _1 recorded at centre-of-mass energies from 88 GeV to 202 GeV, 
a search for two photon decays of Higgs bosons produced in association with a fermion pair 
has been performed in the mass range from up to 120 GeV/ c 2 . No evidence for resonant 
production of photon pairs has been found. In the framework of Higgs bosons with standard 
model coupling to gauge bosons, a 95% confidence level upper limit on the Higgs boson 
branching ratio to two photons has been obtained for Higgs boson masses from to 109 GeV/ c 2 . 
In the fermiophobic model in which the Higgs boson couples exclusively to gauge bosons, a 95% 
confidence level lower limit on the Higgs boson mass has been set at 100.7 GeV/ c 2 . 
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